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Abstract 

The rotor-router model is a deterministic analogue of random walk 
invented by Jim Propp. It can be used to define a deterministic aggre¬ 
gation model analogous to internal diffusion limited aggregation. We 
prove an isoperimetric inequality for the exit time of simple random 
walk from a finite region in and use this to prove that the shape of 
the rotor-router aggregation model in Z'^, suitably rescaled, converges 
to a Euclidean ball in 


1 Introduction 

Given a finite region A C let A be the (random) region obtained by 
starting a random walk at the origin, stopping the walk when it hrst exits 
A, and adjoining the endpoint of the walk to A. Internal diffusion limited 
aggregation (“internal DLA”) is the growth model obtained by iterating this 
procedure starting from the set containing only the origin: Ai = {o}, = 

{An-i)'. Lawler et ah [5] showed that the region rescaled by a factor of 
converges with probability one to a Euclidean ball in as n —>■ oo. 
Lawler nm estimated the rate of convergence. 

Jim Propp has proposed the following deterministic analogue of internal 
DLA in two dimensions. At each site x G A is a “rotor” pointing North, 
East, South or West. A particle is placed at the origin and performs rotor- 
router walkvMiil it exits the region A. during each time step, the rotor at the 
particle’s current location is rotated clockwise by 90 degrees, and the particle 
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takes a step in the direction of the newly rotated rotor. The intent of this 
rule is to simulate the hrst-order properties of random walk by forcing each 
site to route approximately equal numbers of particles to each of the four 
neighboring sites. When the particle reaches a point not in A, that point 
is adjoined to the region and the procedure is iterated to obtain a sequence 
of regions An- For example, if all rotors are initially pointing north, the 
sequence will be begin Ai = {o}, A 2 = {o, (1,0)}, A^ = {o, (1, 0), (0, —1)}, 
etc. 

There has been considerable recent interest in obtaining a shape theorem 
for the rotor-router model analogous to that for internal DLA izma. Much 
of this interest has been driven by simulations in two dimensions, which 
indicate that the regions An are extraordinarily close to circular (Figure Q). 
Despite the impressive evidence for circularity, very little progress has been 
made until now in the way of rigorous results. In one dimension, with rotors 
alternately pointing left and right, the dynamics of the model are simple 
enough to analyze explicitly; in this case the hrst author has shown ^2] that 
the deviation from a ball (symmetric interval) is bounded independent of 
n. In addition, various modihcations and extensions of the one-dimensional 
model are amenable to explicit analysis, and analogous shape theorems are 
known in some of these cases izma. In two dimensions, the hrst author has 
shown ^21 that the region A^, contains a disc of radius proportional to 
In higher dimensions, the model can be dehned analogously by repeatedly 
cycling the rotors through an ordering of the 2d cardinal directions in 
until now nothing was known about the shape for d > 3. 

Denote by i? A S' the symmetric difference of sets R and S. For a region 
A C Z*^, we write A^ for the union of closed unit cubes in centered at the 
points of A. We write C for d-dimensional Lebesgue measure. As a special 
case of our main result. Theorem ini we obtain the following. 

Theorem 1.1. Let {An)n>i be rotor-router aggregation in 7A, starting from 
any initial configuration of rotors. Then as n —> 00 

C{n-^/^A° A B)^0, 

where B is the ball of unit volume centered at the origin in M'^. 

In Theorem o we prove this result in the more general setting of arbi¬ 
trary rotor stacks of bounded “discrepancy;” see section 2 for details. We also 
give an explicit bound on the rate of convergence. Here we should emphasize 
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Figure 1: Rotor-router aggregate of 270, 000 particles. 


that much work remains to be done if one hopes to explain the almost perfect 
circularity found in Figure ^ The form of convergence in Theorem o is not 
as strong as the convergence in the shape theorems for internal DLA CD]. 
In particular, Theorem II .11 does not preclude the formation of long tendrils, 
or of “holes” close to the origin, provided that the volume of these features is 
negligible compared to n. So we hope that this represents only the beginning 
of attempts to understand the shape of the model. 

A major component of the proof is an isoperimetric inequality for the 
expected exit time of random walk from a region in Z'^. Because of its 
intrinsic interest and possible utility in other applications, we state it here. 
Given x E A d let e^iA) be the expected time for simple random walk 
started at x to first leave the region A. We denote by the volume of the ball 
of unit radius in R'^. The following result shows that e^iA) is asymptotically 
maximized among all regions of a given size when A is a ball. 

Theorem 1.2. If X E A and |A| = n, then for sufficiently small e > 0 

ex{A) < eo(R„)(l 0(n“")), 

where Bn is the lattice ball of radius {n/cOdY^^^: 

Bn = {yEZ'^ ludWyW^dn}. ( 1 ) 
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Here o denotes the origin in and || ■ || the Euclidean norm. In Theo¬ 
rem ED we give an explicit bound for the exponent in the error term. 

The proof of Theorem 11.11 also relies on a special “abelian” property of 
internal DLA. Beginning with a conhguration of hnitely many particles in 
suppose that at each time step we choose a site occupied by more than one 
particle and move one particle at that site by one step according to simple 
random walk. After a succession of such choices, the process ends when no 
site is occupied by more than one particle. Diaconis and Fulton ji] discovered 
the remarkable fact that both the distribution of the shape of the resulting 
cluster of occupied sites and the distribution of the number of time steps 
taken to reach the hnal state are independent of the choices made along the 
way. Lawler et ah jU] (sec. 6) found a particularly simple proof by labeling 
the particles and moving only the highest-labeled particle at a given site. 

2 Convergence to a Ball 

For X, 1 / G Z'^ we write x ^ y ii \ \x — y\\ = 1. By a “region” A C Z'^ we 
will always mean a hnite region. We write |A| for the cardinality of A. The 
boundary of A is the region 

dA = {x e I X ~ 2 / for some y G A}. 

Simple random walk in WA will be denoted by Xq, Xi,...; the probability and 
expectation operators for walk started at Xq = x will be denoted and E^,. 
Given a set S' C Z*^ we write 

Ts = mf{j > 0 I X, G 

for the hrst hitting time of S. If A is a region in Z"^ and x G A, we denote 
by ex{A) the expected time for simple random walk started at x to exit the 
set A: 

ex{A) = ExTqa- 

Note that ex{A) = 0 for x ^ A. We write 

e(A) = supe^(A). (2) 

X 

We study the following mild generalization of the rotor-router model in 
Z'^. Fix a positive constant D, the discrepancy. At each site x G Z*^ is 
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an infinite stack of rotors ri,r 2 ,... each pointing in one of the 2d cardinal 
directions. On the i-th visit to the site x, the particle travels in direction rj. 
We reqnire that for any direction £ and any positive integer m, 


< m \ Vi = e} 


m 


< D. 


(3) 


Observe that the original rotor-ronter model with cyclically repeating rotors 
satishes this condition with discrepancy D = 1. 

We will phrase onr resnlt in terms of the qnantity 


0(n) = snp {e{A) - eo{Bn)). 

\A\=n 


Let 

n 

In Theorem lO we show that (j){n) = for snfficiently small e. 

Write C for Lebesgue measnre in For ^4 C we write 


= A + 


1 1 
2’2 


C 


For RcW^ and A G M, we write 


XR = {Ax I X G R}. 


Onr main resnlt is the following. 

Theorem 2.1. Let {An)n>i be rotor-router aggregation in lA using any con¬ 
figuration of rotor stacks with discrepancy at most D. Then 

A B)= ^ q 

as n ^ oo, where B is the ball of unit volume centered at the origin in 
and C = C{d) is a constant independent of n and D. 

Remark. In two dimensions, Theorem Id. II implies that $(n) = hence 
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For d > 3, Theorem ED gives 


$(n) = 

hence ^ 

C{n-^/^A° AB) = 0(n“2^^) + 

Given a hnite region A C dehne 

V’(^) = ikip. 

x^A 

Among regions A (Z with |A| = |-B„|, the qnantity 'ip{A) is minimized 
when A = Bn- The idea of the proof of Theorem 12.II is to show that 'ip{An) 
cannot be mnch larger than 'ip{Bn) iProposition I2.2|) . To estimate 'ijj{Bn), 
notice that 


Biin/udy/^ -d)cB°Z + d), (5) 

where B{r) is the ball of radius r centered at the origin in R'^. From this we 
obtain 

ABn) = ( 6 ) 

Proposition 2.2. 'il^{An) = + *h(n) + _ 

To prove Proposition 12.21 we hrst relate the quantity 'ip{An) to the total 
number of steps T„ taken by the rotor-router walks of the hrst n particles. 
We will make use of the identity 

+ (Xi + 1)2) = 1. (7) 

i=l 

where A^f denotes the discrete Laplacian of the function / at the point x: 

^ E ■'■(») - = iE./(n) - /(i). 

y'^x 

Lemma 2.3. ‘ip{An) < Tn + 8y/dD '^xeAn 11^11 "^dDn. 
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Proof. Given a finite set of particles at locations a:(l),... ,x{n) G define 
the quadratic weight of the configuration to be 

n 

w = w{x{l).,, x{n)) = I |a;(f)| 1^. 

i=l 

At any given time during rotor-router aggregation, each site x has routed an 
equal number of particles to each of its neighbors, plus an error of at most 
2D extra routings to each neighbor y x. By (ED it follows that the net 
effect of the first m routings from a site x is to increase the total weight by 
m plus an error of at most 2D Ill^lP-lblPl- 

Starting with n particles at the origin, let the particles perform rotor- 
router aggregation one at a time until all of is occupied. This process 
involves a total of routings, so the net change in weight is T„ plus an error 
of at most 


2D Y. lllilP-IMI" 

X'^y&An 


d 

= 2D 1| -f \2xi 

X£An 2=1 
d 

< 2D^^(4|i.| + 2). 

XSAn i=l 

< sVdD E I |a;| I -|- AdDn, 

X&An 


1 |) 


where the last inequality is Cauchy-Schwarz. The result now follows from the 
fact that the initial conhguration has weight zero and the hnal conhguration 
has weight fj^An). □ 


We next relate the quantity Tn to the expected exit time ex{An). For any 
region A we have the Laplacian identity 


Axex(A) = -1, X e A, (8) 

We will need an estimate for the exit time eo{Bn) from the discrete ball Bn. 
By stopping the bounded martingale ||Wt|p — t at the time T = Tqb„ when 
the walk exits Bn, we obtain 

eoiBn) = EoT = EoIIXtII' = (9) 


This asymptotic result for random walk becomes exact for Brownian motion 
when Bn is replaced by a ball of radius {n/cOdY^^ in Mf. 
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Lemma 2.4. If |A| = n, then Y.x^y&A\ex{A) - ey{A)\ = 
Proof. By Cauchy-Schwarz, 




< 4n ^ (e^(A) 

X'^yGA 


eyiA)f. 


To bound the latter sum, the fact that ex{A) = 0 for x G dA, together with 
(jHl) and Theorem Id.II implies 


5^ (eAA) - e,(A)r 

xr-^y^AUdA 


2 ^ ex{A){ex{A) - ey{A)) 

X'--'yGA\JdA 

2 ^ e^{A) ^(ej;(A) - ey{A)) 

xGA j/~x 

2'^ex{A){-2dA^ex{A)) 

xGA 

Ad^e^{A) 

xGA 


Lemma 2.5. T„ = neo{An) - Y1 ,x&a^ ex{An) + 

Proof. Given a hnite set of particles at locations x(l),..., x{n) G A^ dehne 
the exit weight of the conhguration to be 

n 

r] = r]{x{l),...,x{n)) = ^e^(j)(A„). 

j=i 

By (0 and 0, the net effect of the hrst m routings from a site x is to 
decrease the total weight rj by m, plus an error of at most 2D Yhyr^^x \^x{An) — 
ey{An)\. Beginning with n particles at the origin and ending when A^ is 
completely occupied, the total decrease in weight thus differs from by at 
most 2D J2x'~~^y&An \^x{.An) — ey{An)\. The result now follows from LemmaEH 

□ 


Lemma 2.6. T„ < _|_ $( 77 ,) -|- OifDrP^^^Af 


Proof. It is here that we will use the abelian property of internal DLA men¬ 
tioned in the introduction. Beginning with n particles at the origin, let each 
particle pj in turn perform simple random walk until it either exits the region 
An or visits a site not visited by any of pi,..., Pj-i- The expected time taken 
by this procedure is at least neo{An) — YhxeAn ^x{An). Letting the particles 
which were stopped upon exiting An continue walking until they also reach 
unoccupied sites, the abelian property implies that 

neo{An) - ^ e^(An) < 

X&An 

where is the expected number of steps taken by the hrst n particles in 

internal DLA. By © we have 

n 

n 

< + 4(n) 

The result now follows from Lemma Ym □ 

We can now prove Propositiou I2.2l bv means of a bootstrapping argument. 
Lemma 2.7. If ^p{An) = 0{D°‘n^) for some a > 1, /3 > 1 + |, then 

^/J{An) = 

Proof. By Cauchy-Schwarz 

[ ^ Ikll I ^ ^ IkiP = (10) 

\x&An / XGAu 

Lemma o now shows that 

f’{An) <Tn + 0(Di+“/2ri(i+/3)/2)^ 

and the result follows from Lemma f2.01 and Q. □ 
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Proof of Pror)osition \2/A Since is connected, ||a:|| < n for all x G 
hence 'ip{An) = 0{n^). The sequences defined by 


aQ — 1, ttm+l — 1 + 


CXy, 


Po 3 , Pm-\-l 


1 + Pr} 


have the explicit forms 


a^ = 2- 2-”^, = 1 + 2 


1—m. 




hence 


1 2 

[log d/log2] < 2 — —, p[log d/log 2] < 1 + ^, 


where [x] denotes the least integer > x. By iteratively applying Lemma f2.71 
we obtain after [log d/log 2] iterations 

p{An) = 


Equation (uni) now gives 


Ikll = 0{D 




L+l/c/^ 


XSAn 

hence by Lemmas 12 .dl and L 

PiA^) < T„ + 0(D2-i/4V+i/d) 


d 


d + 2 

The result now follows from 


+ <h(n) + 




□ 


For the proof of Theorem 12.11 it will be useful to rephrase equation © in 
terms of the radius of the ball: 


( 11 ) 


We will also need a simple estimate for the cardinality of Bn- Recall that 

-d)cB°C B{{n/u;,Y/<^ + d), (12) 


10 





where B{r) is the ball of radius r centered at the origin in It follows that 

\Bn\=n + 0{n^-^/^). (13) 

Expressed in terms of the radius, this becomes 

\B^^^d\ = u)dr'^ + 0{r^~^). (14) 

Better estimates exist (see e.g. 0), but we will not need them. 

Proof of Theorem A2. 1\ We will show that 

\An A ^ 

By (HI this implies 

AB) < C{n-^/'^A° A A B) 

< -\A^AB^\+Oin-^/^) 

n 

which gives the theorem. 

To prove m, we hrst observe that if A i?„| = E, then by (|T!^ we 
have 'f{An) > 'fiA), where A = (B^ \ S^) U 5*+ is the region formed by 
deleting from B^ an outer spherical shell S'- of cardinality E/2 + 
and adjoining an adjacent spherical shell 5'+ of cardinality E/2 — 0(n^“^/'^). 
The outer radius of S_ and inner radius of 5'+ are both equal to r = {n/ojdY^^- 
Solving for the inner radius r_ of S_ and the outer radius r_|_ of S_|_ in terms 
of E, we obtain from 03 

r, = (!i±lVi±^(AAy'Vo(i). 

Since jg ^ convex function of t, it follows that 

r^+2 = Coin ± E/2)1+2/'^ + 0{n^+^/^). (16) 

Equation m yields 

AiAn) - AiBn) > A {{B^^rl \ ^0.^^) U - AiBn) 

= Ci[r++^ - 2r'^+2 + rl+^] + 0(n^+^/'^). 
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Applying (ITT)|l and expanding (n± ± using 

the binomial theorem, all terms involving an odd power of V cancel, and all 
terms involving an even power of V are nonnegative, hence 

Solving for V and applying Proposition 12.21 we obtain 

which yields (fTH|l . □ 

3 Isoperimetric Inequality for Exit Times 

In this section we show that among all regions A C of a given size, the 
maximal expected exit time 


e(A) = supea;(A) 

xGA 


is asymptotically maximized when A is a ball. 

Theorem 3.1. If A dlA and |A| = n, then 

e(A) < eo{Bn) + log^n) (17) 


where 


Id 



d=l 
d = 2 
d>3. 


(18) 


Recall from (0 that eo{Bn) = The case d = 1 of Theorem 13.11 is 

an elementary gambler’s ruin calculation; for the remainder of this section we 
assume d > 2. We have not attempted to optimize the exponent in the error 
term for d > 3, preferring instead to give the cleanest possible arguments. A 
careful optimization of Lemma 13.31 should yield a somewhat smaller error. 

Isoperimetric inequalities of this type have long been known for the exit 
time of Brownian motion from regions in The hrst such result goes back 
to Polya na. who showed that the disc in maximizes “torsional rigidity” 
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among all simply connected plane domains of a given area. Aizenman and 
Simon pQ use a rearrangement inequality of Brascamp et al. jS] to prove 
that a Euclidean ball in simultaneously maximizes all moments of the 
Brownian exit time among all regions of a given volume. 

The proof of Theorem Id. II will proceed in several steps. We hrst appeal 
to a rearrangement inequality of Pruss m to reduce to the case when A 
has a certain weak convexity property fLemma ld.2|l . This convexity enables 
us to estimate the exit time from points close to the boundary by bounding 
the hitting time of an orthant in Z'^ 1 Lemma Id.dj) . The Einmahl extension 
j2j of the Komlos-Major-Tusnady strong approximation [S] (see also |15p 
yields a close coupling of random walk in A and Brownian motion in the 
corresponding region A^ C so that the random walk is likely to be close 
to the boundary of A when the Brownian motion exits A^. Finally, the 
theorem of Aizenman and Simon is used to bound the expected exit time of 
Brownian motion from A^. 

Denote hj ei,... ,ed the standard basis for and by Hi the hyperplane 
spanned by £i,..., £j_i, Sj+i,..., Given a region A C Z'^, for each x E 
let 


ai{x) = #{j eZ \ x + jSiE A}. 


The Steiner symmetrization of A with respect to the hyperplane Hi is the 
region 



X + jSi 


ai{x) ^ ai{x) 1 
2 ^ - 2 ]' 


In words, UjA is obtained by compressing to an interval each column of points 
in A lying above a point x E Hi, and then centering that interval about the 
hyperplane Hi, with preference for the positive side of the hyperplane if the 
interval has even length. In particular, \aiA\ = |A|. 

We say that a region A C Z'^ is orthoconvex ii x E A and x + ksi E A, 
k > imply x + jSi E A for all 0 < j < A;; equivalently, any line in JA parallel 
to one of the coordinate axes meets A in an interval (possibly empty). 


Lemma 3.2. For each n > 1 there exists an orthoconvex region A (Z JA 
which maximizes the quantity e{A) among all regions in TA of size n. 


Proof. Denote by A the set of all connected regions A C Z'^ of size n con¬ 
taining the origin. Clearly, the maximum value of e(A) among all regions 
of volume n is attained by a region A E A. It follows from the difference 
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equation (jS)) and a rearrangement inequality of Pruss m that e does not 
decrease under Steiner symmetrization. On the other hand, the quantity 

d 

xGA i=l 

strictly decreases under Steiner symmetrization unless A is already symmet¬ 
ric. Choosing from among those regions in A which maximize e one which 
minimizes we obtain a region that is Steiner symmetric about every coor¬ 
dinate axis, hence orthoconvex. □ 


If A is orthoconvex, then any point x G dA has a “supporting orthant” Q 
based at x lying entirely outside A. To bound the time to exit A from points 
near the boundary, it suffices to bound the time to hit this orthant. We write 
C{x, r) for the ball of radius r (cube of side length 2r -|- 1) centered at x. 
Simple gambler’s ruin considerations imply that 


ExTdc{x,r) = O(r^). 


(19) 


Lemma 3.3. Let 

1,..., d], and let 


Then if r > 3k, 


Q be the nonnegative orthant {x G Z'’ 


p{k,r)= sup ^x{Tq > TQc(o,r))- 


||a:||oo<fc 


)-d 


p{k,r) < ( 1 - — ) p{3k,r). 


Xi > 0, i = 


Proof. Given x & lA with ||a;||oo < k, let Q{x) = {y & lA \ Pi > Xi, i = 
l,...,d} be the orthant parallel to Q based at x. Subdividing the cube 
C = C{x, 2k — 1) into 2'’* cubes of side length 2k — 1, the intersection Q{x)r\C 
consists of one of the cubes in the subdivision. By symmetry. 


Fx{XT,a e Q{x)) > 2-A 


Now if y is any point in dC fl Q{x), then an entire boundary face of the cube 
C = C{y, k — 1) lies in Q (Figure 2), hence by symmetry 

The result now follows from the observation that the norm of any point 
on the boundary of C or C" is at most 3k. □ 
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Figure 2: Diagram for the proof of Lemma fH.HI 


For a domain D C denote by Tq^ the time when Brownian motion 
exits D. The theorem of Aizenman and Simon PP implies that is 

maximized among domains D of volume n when D is a ball and x is its 
center. Since vol(A'^) = |A| = n we obtain lE^Tg^ < for all 

X G A^. (We adopt the notational shorthand dAP := c}(A'^).) Chebyshev’s 
inequality gives 

sup > 2{n/uJdf/'^) < 

xGA^ ^ 


hence 

> 2m{nlu^dfl^) < 2 -”^. ( 20 ) 

The following is a rehnement of Lemma fH.HI in dimension two. 


Lemma 3.4. In dimension d = 2, 

p{k, r) < c 


there are eonstants a and c such that 

/ k + a\ogr\^^^ 

\r — a logr / 


Proof. Applying the map z i—>■ the conformal invariance of harmonic 

measure for planar Brownian motion implies m that 


PBM{k,r) 


sup 

lhl|oo<fc 


P frriBM 

dC{o,rp 




( 21 ) 
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By the strong approximation theorem isiEiiini there exists a constant a > 0 
and a coupling of simple random walk in l/' with Brownian motion in W^, 
so that, except for an event Ei of probability at most the coupled paths 
are separated by a distance of at most alogr — 2 up to time s = logr]. 
Let E 2 be the event that the Brownian motion has not exited C{o,r)^ by 
time s. By fOUj) we have Fx{E 2 ) = 0{l/r). On the event E^ 0 if the 
random walk exits C{o,r) before hitting Q, the Brownian motion must exit 
C{o,r — a log r) before hitting the translated quadrant + {a log r, a log r); 
hence 


p{k,r) < pBuik + alogr,r - alogr) + F^{Ei) + F,^{E 2 ). 

The result now follows from (j2^- □ 

Proof of Theorem, \S . /L Denote by i ?3 the event that the random walk and 
Brownian motion paths in the strong approximation coupling are separated 
by distance more than 61ogn — d before time s = \rP^^logn]. Choosing b 
sufficiently large we can take F^E^) < K Write T = and denote by E^ 

the event that T > s. By fOUll we have Fx{E 4 f = 0{^) for all x G A^. 

On the event E^ fl E^ the location Xt of the random walk when the 
Brownian motion exits A is distance at most b log n from dA. Let Q (Z A'^ he 
the supporting orthant at a point Y G dA within distance 61ogn of Xt- For 
j > 1 let Ej be the event that after time T the walk travels to distance 
3-^6 log n away from Y before hitting Q. Iteratively applying Lemma l!01 with 
initial value k = b log n we obtain 

Write m = \ • By (j221) we have 

F^{E^) = 0{n-'^'^), d>3. (23) 

In dimension two, Lemma 13.41 with k = 61ogn and r = 3-^61ogn gives for 
j <m 

Px(F;) < (24) 

Taking j = m we obtain Fx{Em) = Thus (12 3|) holds in dimension 

two as well. 
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On the event E^nElT} Ej, the time for the random walk to exit A is at 
most the exit time for Brownian motion plus the time for the walk to go an 
additional L°° distance (3-^ + 1)6 log n. Hence 

m 

TdA <T + ^ lFj^9C(XT,(3J+l)felogn) + + l£;3(s + f^) + 1 ^ 4(5 + Ta) , 

i=i 


where T is the additional time taken to exit A if the walk travels to distance 
3"*61og?7, from Y before hitting Q] and T) for i = 3,4 is the additional time 
taken to exit A after time s if the event Ei occurs. Taking expectations and 
applying (UHl), we obtain by the strong Markov property 


Ex Tqa < 


\ 2/d m 

— + Cl V ¥^{Ej)3^^ log' n + E„ T + — -f25) 

J 


i=i 


2s 


n 


O 


n 


(E^Ts + Ea,T4). 


By (j22I), (1^ and (OH) . 


5^Px(C,) 3'^' < CaP. (F^)3'”^ = 0(n'/"-^Mog'n). 
i=i 


Maximizing (123) over all a; G H we obtain 


e{A) < 


J 

and solving for e(y4) yields 


+ 0(n'/'^ log'n)+ 0(n '^'^)e(H)+ 0(n'/'^ ^logn) + —e(H) 

n 


/ \ 2/d 

e(H) < f — j + log'n). □ 
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